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The elementary excitations in antiferromagnets are magnons, quasiparticles 
with integer spin and Bose statistics. In an experiment their density is controlled 
efficiently by an applied magnetic field and can be made finite to cause the 
formation of a Bose-Einstein condensate (BEC). Studies of magnon condensation 
in a growing number of magnetic materials provide a unique window into an 
exciting world of quantum phase transitions (QPT) and exotic quantum states. 

INTRODUCTION 
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Not long after Bose and Einstein described quantum statistics of photons 
Bloch apphed the same concept to excitations in sohds. He treated misahgned spins in a ferro- 
magnet as magnons, particles with integer spin and bosonic statistics, to describe the reduction 
of spontaneous magnetisation by thermal fluctuations [3]. Matsubara and Matsudapointed out 
an exact correspondence between a quantum antiferromagnet and a lattice Bose gas [4]. It is thus 
natural to ask: can these bosons undergo Bose-Einstein condensation and become superfluid? The 
answer is yes. The question of order in spin systems and possibility of BEC of magnons has been 
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investigated theoretically for a number of quantum antiferromagnets 
Although in simple spin systems the concept of BEC can be applied, in practice factors such as the 
large value of exchange constants and the presence of anisotropies violating rotational symmetry 
may restrict their usefulness. However, the analogy between spins and bosons has shown to be 
very fruitful in those antiferromagnets where closely spaced pairs of spins 5 = 1/2 form dimers 
with a spin-singlet (5 = 0) ground state and triplet {S = 1) bosonic excitations called variously 
magnons or triplons. In such systems BEC has been predicted to occur [7] and was experimentally 



observed [8 



14 ] in the dimer system TlCuCls. This started a flourish of activity on the subject 



and the hunt for such transitions in other magnetic materials. 

Other examples of condensates from elementary excitations were recently found in magnetic thin 
films lisj] and semiconductor microcavities [3]. Lattice bosons and several transitions occurring in 



these systems have been studied extensively on uhra-cold atomic gases in optical lattices [l^, Q, 
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20( 1 . Here we present an overview of recent developments in the field, and discuss how quantum 



antiferromagnets offer several advantages that set them apart from the other model systems in 
which the phenomenon of BEC occurs. We also discuss ways to go beyond the physics of simple 
BEC condensation and to look at the fascinating new quantum phases of interacting bosons on a 
lattice. 

BOSONS IN MAGNETS 



Let us illustrate the basics of i nag non BEC in real dimerized antiferroinagnets, such as ex- 
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291] and BaCuSi206 



36( 1 . The lattice of magnetic ions in such materials (Fig. la) can be 
visualized as a set of dimers, pairs of copper ions Cu^^ carrying 5 = 1/2 each and interacting via 



Heisenberg exchange: 



H — ^ Jo Si^i ■ 82,4 + ^ Jmnij 

i {mnij) (ni) 



(1) 



where H denotes an external magnetic field in z-direction, i,j number dimers, and, m, n = 1,2 
their magnetic sites. The intra-dimer exchange is the strongest interaction, which happens to be 
antiferromagnetic, Jq > 0, so that an isolated dimer has a ground state with total spin S = and 
a triply degenerate excited state of energy Jq and spin S = 1 (Fig. Ic). It is convenient to identify 
the triplet state with the presence of a triplon, a bosonic particle with S = 1, and the singlet state 
with the absence of a triplon, see Fig. lb. As long as mter-dimer interactions are weak, the ground 
state consists of non-magnetic singlets. It is disordered down to absolute zero temperature without 
long-range magnetic order. The triplon excitations are made mobile by effective and weak mter- 
dimer couplings J1.2,... > 0, from a sum over single-ion interactions Jmnij, see Fig. lb. For dimers 
forming a square lattice (for simplicity) the energy of a triplon with spin projection 5^ = 0, ±1 is 



e(k) = Jo + Ji[cos (kxtt) + cos (kya)] — qhbHS^ , 



(2) 



where k = (kx, ky) is the particle wavevector, a is the lattice constant, and D = 4Ji the bandwidth, 
see Fig. Ic. The energy-momentum dependence of the triplons and singlet-triplet correlations have 
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been measured directly by inelastic neutron scattering 

The Zeeman term —gfisHS^ controls the density of triplons. As the magnetic field increases, 
the excitation energy of triplons with = +1 is lowered and eventually crosses zero, as shown in 
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Fig. Ic and 2a. This defines two critical magnetic fields Hd and Hc2 in the phase diagram, see 
Fig. Id. At zero temperature, below Hd the magnetisation mz{H) is zero and only singlets exist. 
Between Hd and Hc2 the magnetisation increases as the triplon band fills up, see Fig. Ic. Above 
Hc2 each site is occupied by a triplon and the magnetisation saturates at one per dimer. 

The bosonic nature of triplons is guaranteed by the simple fact that spin operators of two 
different dimers commute. However, because a dimer can hold at most one triplon, the bosonic 
picture requires the introduction of a hard-core constraint to exclude states with more than one 
quasi-particle per dimer. The constraint, which can be interpreted as a strong short-range repulsion 
between the bosons, poses a difficult theoretical problem. But close to Hd their density is small 
and collisions between bosons are rare; in this limit interactions can be fully taken into account. 
By particle-hole symmetry, a similar simplification occurs near Hc2- 



BEC OF TRIPLONS 



The nature of the ground state above Hd, its interpretation as a BEC of magnetic quasi- 
particles, and the tuning of the particle density become particularly accessible if the spin Hamil- 
tonian in Equ. ([T|) is rewritten in the second-quantized form 0, ^ 

H = Y.^'^o - h)a\ai + ^ Uja\aj + \Y^ Uija\a]ajai (3) 



where a\ (aj) create (annihilate) a boson on dimer z, and h = gfisH is the effective field. The 
term Uj defined by the transverse component of the inter-dimer interactions Jf'^ etc describes 
hopping between sites i and j thus endowing the triplons with kinetic energy, see Fig. lb. Uij 
from the longitudinal components Jf etc is the repulsion energy arising when two triplons occupy 
neighboring sites i and j. The density of triplons is directly controlled by the magnetic field that 
acts thus as a chemical potential for the bosons [38 1. 



Below Hd the ground state is a quantum-disordered paramagnet formed by the singlet sea 
(triplon vacuum) and can be approximated by the direct product of singlet states on each dimer, 
= IS, S^)i with S = = 0. Once the spin gap is closed at Hd, a Bose condensate is formed. 
Since the bottom of the triplon band is located at a nonzero wavevector ko = (^, ^) (Fig. Ic), 
the wavefunction of the condensate varies in space as exp (iko • r) . In this phase the state of an 
individual dimer is well approximated by a coherent superposition of the singlet and the = +1 
triplet: 

field [7|,lll|,l]J. 



= ai(H)\0, 0)j + Pi{H)\l, where the amplitudes and Pi depend on the magnetic 
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In the spin language, the condensate corresponds to a spontaneous magnetic order formed by 
the transverse spin components (Sf) and {Sf), violating the remaining rotational 0(2) symmetry 
of the Hamiltonian ([T]). To make the analogy with the traditional BEC manifest, one can form a 
U{1) order parameter (Sf+iSf). The phase corresponding to the angle of the spin in the XY plane 
is thus the phase of the wavefunction in the boson language. At Hd the paramagnetic phase at 
low fields makes a transition into a canted antiferromagnet with long-range magnetic order in the 
plane perpendicular to the field (Figs. Id and 3a). The staggering of the transverse components 
of magnetisation reflects a nonzero wavevector of the condensate ko. The critical properties of the 
magnet in the vicinity of this phase transition are governed by the quantum critical point (QCP) 
of the BEC universality class located at T = and H = Hd- Indeed close to Hd the bosons are 
extremely diluted, and the effects of the interaction become weak, despite its hard-core nature. 
Close to the QCP the phase boundary Tc{H) follows Q] a power law Tc oc {H — Hd)'^ with a 
universal critical exponent cf) = z/d, which depends only on the dimensionality (d) and dynamical 
critical exponent z = 2 for a quadratic triplon energy band. The upper critical dimension for the 
QCP is dc = 2. 

The BEC QPT has been observed in a growing number of dimer based magnetic insulators, such 



as ACuCla (A=T1, K, NH 



and (C4Hi2N2)Cu2Cl6 



BaCuSiaOfi, Cu(N03}2i2.5D20, Cs3Cr2Br9, (CH3)2CHNH3CuCl3, 
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47( 1 . The physics of triplon 



condensation presented for the dimer system remains essentially unchanged since these QCPs are 
in the same universality class. Quasi-one-dimensional arrangements of S* = 1 moments in spin 
chains, e.g. Haldane chains, and even field-saturated frustrated antiferromagnets can be described 
within the same framework, as was successfully done for nickel-based materials and CS2CUCI4, 



respectively 
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Experimentally, the static and dynamic magnetic properties have been studied at and around 
the BEC QCPs in such materials by many experimental techniques. The evolution of magnetic 
excitations across the QCP at ff^i in a three-dimensional network of dimers was studied in TICUCI3, 



see Fig. la, 2a, and 2c 
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22l |. The softening of the triplet mode = +1 in the quantum 



disordered phase. Fig. 2a, is followed by a dramatic change in the nature of the excitation spectrum 
above the critical field Hd- In particular, as can be expected in a system with a spontaneously 
broken XY symmetry in the plane perpendicular to the applied magnetic field, a 'Goldstone' mode, 
with a dispersion linear with momentum appears 
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Historically, the temperature dependence of the magnetisation mz{T) at fixed field H > Hd 
in TICUCI3 (Fig. 2d) and the scaling of the critical temperature Tc{H) were the key experimental 



observations that confirmed the theoretical model [3, of the magnon BEC in dimerized quantum 
antiferromagnets. The minimum and cusp in mz{T) at the finite-temperature transition can not 
)e explained within simple mean- field theory jsS], but are a consequence of magnon condensation 
3, B] and as such also occur in lower-dimensional magnets. For example comprehensive data from 
bulk experimental measurements establish the phase diagram and critical exponents for BaCuSi206 
(Fig. 2b). The temperature crossover exponent is found to be (/> = 2/3, as expected for a BEC QCP 
in three dimensions. However, below 1 K and down to 35 mK the phase boundary unexpectedly 
becomes linear, indicating that < 
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36( 1 . We discuss this anomaly below. Even in the 



case where a full condensate does not exist, such as in one-dimensional systems, similar minima in 
the magnetisation can occur, but now as a simple crossover 541 ]. 



COMPARISON WITH OTHER BOSON SYSTEMS 



The correspondence between a Bose gas and an antiferromagnet is summarized in Table [H 
Despite many similarities, there are also a few important differences between triplons in a magnet 
and atoms in a Bose gas 
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m M, lia, i2o(. mber of atoms is usually 

controlled directly (microcanonical ensemble), while the number of triplons is typically set by the 
magnetic field acting like a chemical potential (canonical ensemble). 

Considerable differences also exist at the practical level. Triplons are much lighter and have a 
much higher density than atomic gases. As a result, condensates survive to much higher temper- 
atures: kelvins in magnets, even room temperature in some cases jl^, as opposed to nanokelvins 
in atomic gases. As solid-state systems, they also allow for a variety of static and dynamic probes 
(magnetisation, specific heat, NMR, neutrons, etc.). On the other hand, the cold atomic sys- 
tems allow potentially to study some out-of-equilibrium situations that are hard to maintain in a 
solid-state device since the condensate is still strongly coupled to the dissipative environment. 

From the point of view of realization, the cold atomic systems have a high degree of control 
and tunability in terms of the interactions and structure. Indeed the structure of the lattice and 
the hopping of the bosons can directly be controlled by varying the strength of an optical lattice. 
The interaction can also be controlled via a Feshbach resonance, allowing in principle to realize 
one's pet Hamiltonian. In spin systems the parameters can be changed only in a limited way by 
the application of pressure or by changing the chemical composition. Efforts in quantum chemistry 
produced good realizations of one-, two-, and three-dimensional dimer materials, e.g. TlCuCls, 
BaCuSi206, and (C5Hi2N)2CuBr4, respectively. In addition, by the very nature of the mapping 
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from spins to bosons, the spin systems offer a definite advantage in reaching the hmit of strong 
on-site repulsion as well as the effects of interactions between nearest neighbors (see the next 
section). The spin systems are therefore an optimal starting point to study situations for which 
these ingredients are important. 

The cold atomic systems have the important advantage that the phase U{1) symmetry is exact. 
In the magnets, the corresponding 0(2) symmetry in the plane perpendicular to the magnetic 
field can be broken by weak anisotropic interactions (dipolar interactions, crystalline anisotropies, 
spin-orbit coupling). Even when symmetry-breaking terms are weak, they become important at 
low temperatures and modify the physics in the vicinity of the QCP. The physics of BEC is also 
altered in the presence of coupling between spins and lattice distortions. Fortunately, it is possible 
to use our good control and knowledge of the pure BEC system to treat theoretically these weak 
deviations as well [2a, l28l . l55l . 1561 ] . Finally the spin systems are excellent to study the critical 
behavior around the quantum critical points, since the boson density can be finely controlled by 
the external magnetic field, and the system is fully homogeneous. Intrinsic density inhomogeneities 
induced by the confining trap in cold atomic systems, make a similar study more difficult. 



TABLE I: Correspondence between a Bose gas and a quantum antiferromagnet. 



Bose gas 


Antiferromagnet 


Particles 


Spin excitations (S = I quasi-particles) 


Boson number N 


Spin component 


Charge conservation U(l) 


Rotational invariance 0(2) 


Condensate wavefunction {ip{r)) 


Transverse magnetic order {Sf + iSf) 


Chemical potential /i 


Magnetic field H 


Superfluid density ps 


Transverse spin stiffness 


Mott insulating state 


magnetisation plateau 



BEYOND SIMPLE BEC 

The interplay of an atomic lattice and strong boson interactions leads to new phenomena beyond 
the simple BEC paradigm. When the boson density approaches a value commensurate with the 
lattice periodicity, one often finds magnetisation plateaus where the triplon density stays constant 
in a finite range of magnetic fields at T = 0. We already discussed two simple examples of such 
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plateaus for H < Hd with zero magnetisation and for H > Hc2 where each site is occupied by one 
triplon (Fig. Id). These simple niz = or 1 plateaus correspond to states with an integer number 
of bosons per unit cell of the lattice. More generally, since in spin systems the S^^S^, terms 
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translate into triplon repulsion, magnetisation plateaus at fractional fillings can also occur 
Such intermediate plateaus arise for strong enough repulsion between spin excitations on adjacent 
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lattice sites and beyond, and correspond in boson language to Mott-insulator phases 
In the simplest case (square lattice, repulsive nearest-neighbor interactions) triplons may form an 
incompressible state by 'crystallizing' in a checkerboard-like pattern (Fig. 3c). Unlike the ground 
state of an integer magnetisation plateau (Fig. 3b), such a state violates the translational symmetry 
of the lattice, and the spontaneous order is characterized by an discrete order parameter {Z2 in the 
case of the square lattice) . Fractional plateaus require strong magnon interactions (in comparison 
to the kinetic energy) and are therefore less robust then integer ones. Nonetheless, plateaus at 
ma gnon fillings 1/3, 1/4, 1/8, and recently for other fractions have been observed in SrCu2(B03)2 
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6d, 



63l |. In this material the geometrical frustration suppresses the kinetic energy of magnons 



thus making magnon repulsion more pronounced [64!]. 

Outside the plateaus, magnons can be freely added to the condensate ensuring a continuously 
varying magnetisation along the applied magnetic field. It has transverse magnetic order violating 
the continuous 0(2) symmetry, is therefore gapless, and corresponds to the super fluid phase of 
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bosons. Its excitations are Goldstone modes, characteristics of this superfluid phase 
incompressible plateau state and the gapless state with a continuously varying magnetisation are 
again separated by a QCP. 

While the dimensionality of an antiferromagnet is determined mostly by quantum chemistry, it 
can also change on the fly as the magnet is cooled down. The "Han purple" BaCuSi206 contains 
layers of dimers with substantial exchange couplings between the layers. Yet the critical behavior 
exhibits a surpri sing crossover from three-dimensional at high temperatures to two-dimensional 



close to the QCP 



33|, 



361 ]. Two effects contribute to this dimensional reduction at low temperatures. 
The interplane triplon hopping is suppressed by geometrical frustration: it vanishes exactly at the 
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wavevector of the condensate 

E, 35, 3 



34( 1 . Additionally, there are inequivalent layers with different 



values of the spin gap 

Another spectacular example of reduced dimensionality is provided by (C5Hi2N)2CuBr4 
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and related organo-metallic compounds. These materials can be thought of as a collection of chains 
of dimers nearly decoupled from each other. As a result, a large part of the H — T phase diagram 



can be understood by starting in the one-dimensional limit [7|, |6 
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71[. At low energies 



the triplons on a single chain form a Luttinger hquid [72]. Such a compound is thus particularly 
suited to address the interesting question of the dimensional crossover between a one- and three- 
dimensional system. In the former, the spins will more behave as fermions, as is well known from 
the Jordan- Wigner transformation, and in the latter, at temperatures and fields for which residual 
exchange coupling between the ladders becomes important, the spins are much more behaving as 
"regular bosons" [73]. The excellent degree of control over the density by the magnetic field as 
well as the various probes providing access to the dynamical spin-spin correlations make quasi- 
one-dimensional magnets ideal for studying the dimensional crossover, a phenomenon of general 
importance for several major areas in physics [t^. 

We would like to mention the tantalizing possibility of finding an exotic phase known as the 
supersolid tJ]. A hybrid of a solid and a superfluid, this phase violates both the translational 
symmetry by forming a density wave and the ?7(1) phase symmetry by exhibiting transverse mag- 



m, 



78|, 



7711 and its existence is under 



791]. Such a phase could oc- 



netic order (Figs. 3b-c). It has been conjectured long ago 
intense debate at the moment following experiments in helium 
cur at the end of a fractional plateau because the transition must accomplish two tasks: create 
a superfluid and destroy a solid (i.e. restore the broken translational symmetry). According to 
Landau's theory of phase transitions, such transitions are either discontinuous or happen in two 
stages. Several theoretical calculations show that supersolids exist in realistic magnetic models 
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The observation of such a phase would thus be an interesting realization 
of the physics of interacting bosons on a lattice. 

Last but not least, advances in growth techniques may allow in the future to study the effects 
601 . 186l . I87f ] and the Bose glass by creating bond disorder, and the influence 



of disorder on bosons 



891 ] . Other directions include the thermodynamics 



of impurities on the quantum-critical states 

of strongly interacting bosons {stI. IooI and quantum coherence at the mesoscopic scale 9ll |. 



CONCLUSIONS AND OUTLOOK 



Quantum spins and quantum dimer systems offer remarkable opportunities to study phenomena 
related to the Bose-Einstein condensation of interacting quantum particles. The high density and 
low mass of spin excitations, magnons or triplons, leads to a robust nature of the condensate, 
which survives to temperatures as high as 10 K. The availability of numerous model magnets and 
physical probes has enabled a detailed study of the critical phenomena and magnetic properties in 
the vicinity of the BEG quantum phase transition. The spin systems have thus proven to be nicely 
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complementary to other systems such as the cold atomic gases for the investigation of BEC. 

In addition to providing a direct connection with the remarkable BEC phenomenon, the boson 
picture offers an intuitive understanding for complex quantum properties of matter that are much 
less physically transparent in the original spin language, and provides access to a cornucopia of 
exciting new problems and questions. For example, a classical description as canted antiferromag- 
nets of the dimerized magnetic materials discussed here in the field-induced order phase misses 
up to 100% of the relevant spin physics close to the intrinsic quantum phase transition. Here 
the possibility to fine tune the density of bosons by using the magnetic field provides a tool to 
determine the role of the interactions in such systems and a wide range of lattice geometries and di- 
mensionalities is available for study. In combination with the large number of experimental probes 
for spins this allows to challenge our knowledge of exotic phases of strongly interacting quantum 
particles, such as BEC in various dimensions, Luttinger-liquid physics, commensurate solids with 
a fractional number of bosons per unit cell, and supersolids combining superfluidity with a broken 
translational symmetry. 
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momentum k„ 
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FIG. 1: BEG of magnons in dimerized quantum antiferromagnets: (a) Dimers in the real material TlCuCla 

with S — \/2 from Cu^+ ions and superexchange via Gi- Q, Q, Q, Q, H Q- 

(b) Dimers on a 

square lattice with dominant antiferromagnetic intradimer interaction Jo and interdimer interactions Ji. 
Triplet states (gray, top) are mapped on quasi-particle bosons (triplons, bottom), (c) Zeeman splitting of 
the triplet modes with gap A and bandwidth D at kn — (-,-). Dispersion of triplons at the critical field 

nnnnn 

Hci |11LI12|.|14|.|21L|22|. (d) Resulting phase diagram with paramagnetic (PM), quantum disordered (QD), 
field-aligned ferromagnetic (FM), and canted-antiferromagnetic (XY-AFM) phase, where BEG of magnons 
occurs. Close to Hd and Hc2 the phase boundary follows a power-law Tc oc {H — Hd)'^ with a universal 
exponent — 2/3 for BEG of magnons 0, Q]. Magnetisation curve mz{H) for 3D dimer spin system with 
plateau at niz = 0, 1 (gapped) 
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FIG. 2: Experimental results on magnon BEG: (a) Zeeman splitting of the triplet modes in TlGuGls up to 
H > Hci measured by inelastic neutron scattering [l4, 27, 3]. (b) Phase diagram of BaGuSi206 measured by 



Mm- 

Dimensional reduction was reported 



torque magnetisation, magneto-caloric effect, and specific heat 
in this material with a crossover from the 3D-BEG critical exponent = 2/3 to = 1 for 2D at temperatures 
close to the QGP [33]. (c) Excitations in the BEG of triplons realized in TlGuGla Q, H, 28|. Goldstone 
mode with linear dispersion around kg. Spin anisotropy generally leads to a spin gap in real materials 
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28 



29|. (d) Temperature-dependence of the magnetisation mz{T) in TlGuGla for fixed magnetic field 
i/, as indicated js]. Minima at the finite-temperature phase transition (vertical arrows), as expected for a 
BEG of triplons. 
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FIG. 3: Lattice boson phases: (a) Canted antiferromagnetic state at Hd < H < Hc2 in a BEC of magnons. 
Longitudinal m^, see Fig. Id, and transverse {XY) magnetisation mxy (b) Triplon condensate in superfluid 
phase. The hopping amphtude t of the quasi-particles leads to a uniform mixture of singlet and triplet 
states, (c) Analog of the Mott-insulating phase for triplons with magnetisation plateau at iriz — 1/3: 
triplons "crystallize" to form a superstructure. This phase is partially observed in SrCu2(B03)2, where 
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quasi-particle hopping is strongly suppressed by geometrical frustration 
a square-lattice of dimers with exchange anisotropy as proposed theoretically 
geometries 
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8C 
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6j|. (d) Supersolid on 



, and for other lattice 



. Both translation and U{1) symmetry are spontaneously broken, i.e. coexistence 



of superfluid and " crystal" of lattice bosons. 



